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Abstract. Recently [A. Xuereb, et at, Phys. Rev. Lett. 105, 013602 (2010)], we 
calculated the radiation field and the optical forces acting on a moving object inside a 
general one-dimensional configuration of immobile optical elements. In this article we 
analyse the forces acting on a semi-transparent mirror in the 'membrane-in-the-middle' 
configuration and compare the results obtained from solving scattering model to those 
from the coupled cavities model that is often used in cavity optomechanical system. 
We highlight the departure of this model from the more exact scattering theory when 
the reflectivity of the moving element drops below about 50%. 
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1. Introduction 

The nontrivial interplay between the external (motional) or internal degrees of freedom 
of a mobile scatterer coupled to a cavity field, and the cavity field itself has attracted 
considerable attention over the past two decades. Use has been made of a cavity 
field to, e.g., interact with single atoms [U El El H], cool atomic motion (5J El IB], 
impose spontaneous order through a Dicke phase transition in an ultracold atomic 
medium [HI [10], couple to the motion of mechanical oscillators [HI H21 [131 HH 115] - 
and even cool this motion down to the vibrational ground state [T6l IT71 [18]. The 
description of these systems, along with most of cavity QED (CQED), follows down the 
path of the 'good cavity' approximation [19]: the cavity mirrors, be they fixed [121 120] 
or moving [2U E] , bound a region of space such that the electromagnetic field in that 
region is cut off from the outside world. An alternative approach, based on a scattering 
picture, is possible. Such an approach can treat very general configurations in one 
dimension, owing to the power of the transfer matrix method (TMM) [221 1231 121] ■ 
In the right limits, the two approaches must of course give rise to the same physics, 
and indeed they do, even in the case of moving boundaries [24]. However, there is no 
guarantee that one TMM model is always equivalent to the same CQED model; it is 
the purpose of this paper to use the specific example of a scatterer inside a cavity, i.e., 
the 'membrane-in-the-middle' scheme [T2| [251 126] to highlight the differences between 
these two approaches. 

Indeed, suppose we have a scatterer, say an atom or a membrane, of reflectivity r 
(0 < |r| < 1) placed inside a cavity which, on its own, can be described very well 
using the 'good cavity approximation.' One of two limiting descriptions is generally 
appropriate for this situation in the CQED picture, (i) If the scatterer were, e.g., an 
atom, with |r| <C 1, the shape of the mode functions of the field inside the cavity will 
not change appreciably. In this case it is valid to treat the atom in a weak-coupling 
approximation and assume that it essentially couples to the unperturbed cavity field||] 
(ii) On the other hand, if the scatterer were a good mirror, with |r| approaching 1, this 
description is no longer valid. Not only does the mirror perturb the shape of the cavity 
field, but in the good-cavity approximation it defines two new modes that communicate 
by tunnelling of photons through the good mirror. This simple example shows the 
power of the TMM approach: the same TMM model is valid for both situations, and 
indeed for any situation in between, including absorbing scatterers, with the value of 
the polaris ability ( of the scatterer determining which of the two situations is being 
described. 

There exists a further, and more fundamental, difference between the two 

| This 'weak coupling' criterion is not related to the so-called strong coupling condition of CQED, 
which refers to the regime in which the internal coherent atom-light coupling leads to a dynamics on a 
time scale shorter than the characteristic decay time of the dissipation processes. This kind of strong 
coupling can be achieved without distorting spatially the empty cavity mode functions of the radiation 
field. 
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Figure 1. The model we consider in this text, drawn schematically. A scatterer S 
interacts with two 'general optical systems' in one dimension, composed of immobile 
linear optical elements, one to either side. S and these two systems are each represented 
by a 2 x 2 matrix. 



approaches. The TMM deals with moving boundary conditions in a way that goes 
beyond merely having a dynamically-changing detuning. Indeed, the mode functions 
themselves in the TMM change dynamically. The implications of such a dynamical 
situation will not be a concern in the following, and we refer the reader to the recent 
work by Cheung and Law |27j for a more thorough discussion of this point. 

The remainder of this paper shall be organised as follows. In the next section we will 
briefly summarise the general solution to the TMM with one moving scatterer (28] [29] . 
The following section will apply this general solution to the study of the 'membrane-in- 
the-middle' model and compare it to the commonly used CQED model [23], following 
which we will conclude. 



2. General solution to the TMM with a moving scatterer 



2.1. Force acting on the moving scatterer 

Consider the generic situation sketched in Fig. [I] Within the TMM, every scatterer 
in the situation is represented by a 2 x 2 matrix M. Free-space propagation at a 
wavenumber k is represented by 
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For a static scatterer, M is related simply to the amplitude reflectivity r and 
transmissivity t of the scatterer, via its polarisability (, which may depend on k: 
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Static scatterers do not change the frequency of transmitted and reflected light. A 
moving scatterer, however, Doppler-shifts reflected light, and we represent this process 
by transforming M into a frequency-dependent operator M f2M [25] . At first order in 
the velocity v of the scatterer, this transformation is remarkably simple and we may 
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write down the general solution for the velocity-dependent force acting on the scatterer 
in closed form [28| [29] . In terms of the notation in Fig. [TJ we can define 



M = Mi x M s x M 2 
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Assuming that the pumping field is monochromatic about some wavenumber ko, B\ = 
Bq 5{k — k ) and C T = C 5{k — k ), we can write the field amplitudes A = f A(k) dk 
and B — J B(k) dk which are given, to first order in v/c, by: 
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where the derivatives are all evaluated at k = ko and act on the frequency-dependent 
terms arising from free-space propagation or a fc-dependent polarisability. 
To obtain these expressions one first solves for A\(k) and D T (k) in terms Bq and Co, and 
then substitutes the results into the matrix equations to obtain explicit expressions for 
A(k) and B(k). Upon noting that these expressions are valid to first order in v/c and 
that the pumping field is monochromatic, the integrals can easily be performed to yield 
Eqs. Q and Q. For single-sided pumping (e.g., Co = 0), these expressions simplify 
significantly. We shall find it useful to express these results in the form A = Ao + \A\ 
and B = Bo + -B\, with Ao,i and Bo,i being independent of v. For conciseness, let us 
now assume that ( does not depend on k. Then, using the elements of Ms, we obtain 



C = j C(k) dk = [(1 - i()A - KBo] + I [(1 - K)A X + 2i(Bo - K&i\ , 



and 
V = 



D(k) dk = [KAo + (1 + K)Bo] + l[2z(Ao - KAt - (1 + z'C)-Bi] 



(8) 



(9) 
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We denote the velocity-independent parts of C and T> by Co and Vq, respectively. The 
force acting on the scatterer can be finally written down as F = F + -Fi, where 
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f (|C| 2 + Im{C})Re{A^}-2Im{C} Re{A^} 
+ (|C| Z - Im{C})Re{B ^} + Im{C} Re{A Bl} 
+ Re{(|C| 2 + ^Re{C})A^}] ; I 
the quantity dF/dv = cF\ will henceforth be called the 'friction coefficient'. 

2.2. Momentum diffusion experienced by the moving scatterer 

The field amplitudes calculated in the previous section related to classical 
electromagnetic fields. We may now impose a canonical quantisation on these fields [24J , 
promoting each field variable A, say, to an operator A, such that (A) = ^2e S/ (h~k ) A, 
S being the mode cross-sectional area. The only two a priori independent modes in 
our system are the two input modes B\ and C T , whose operators obey the usual bosonic 
commutation relations 
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The commutation relations between each of the four fields A, B, C, and D can then be 
built up; because F is correct up to first order in v/c we only need to evaluate expressions 
to zeroth order in this section. The fluctuations in these fields will lead to a diffusion 
in momentum-space, quantified by the diffusion coefficient D. Another contribution 
to D is due to lossy scatterers: any absorptive scatterer effectively couples the system 
to a further, 'loss,' mode that is independent of the input fields and is necessary to 
preserve the canonical commutation relations [21]. Such loss modes can be included 
self-consistently into the TMM |29j. Putting all of this together we can write 
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Knowledge of D and F then allows us to obtain the equilibrium temperature to which 
the scatterer will tend to: 

k B T = -D/{cF x ) , (14) 

where k-Q is Boltzmann's constant. These quantities, which can thus be fully determined 
from our scattering model, are some of the more important quantities of interest 
in optomechanical setups and atom-CQED, and allow us to describe the dynamical 
behaviour of such systems. 
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Figure 2. Our model for the 'membrane-in-the-middle' geometry: the general optical 
systems in Fig. [T] have been replaced by identical mirrors that form a cavity around 
the moving scatterer. We will only consider situations where |x| -C L c . 



3. 'Membrane-in-the-middle' model 

We begin by modelling the system in Ref. [12] : a two-mirror Fabry-Perot cavity with a 
micromirror near its centre, operating at a wavelength A = 1064 nm and having a length 
L c = 6.7 cm, cf. Fig. [2j The micromirror is modelled by its polarisability ( which, in 
light of the small losses observed in practice, is taken to be real and negative. Whereas 
the real experimental system corresponds to |£| < 1, we allow £ to vary freely in our 
model. The two quantities of interest in this section are the intensity of the field close 
to the micromirror, and the friction coefficient acting on the micromirror. The former of 
these gives us knowledge of the resonant frequencies of the cavity and, therefore, of the 
optomechanical coupling, to all orders, between the cavity field and the micromirror. 
The latter is useful in optomechanical cooling experiments; the interest here lies in 
the fact that cooling the motion of a micromirror is one way towards achieving higher 
sensitivity in metrology applications, most notably in gravitational-wave detectors |21j, 
force sensors [30], and magnetometers |31j . 

These quantities are summarised in Fig. |3j with the left panels showing the intensity 
at the mirror and the right panels the friction coefficient acting on the mirror. Each 
subfigure (a)-(f) explores a different value for (. For \(\ <C 1, the cavity field is close to 
the bare-cavity field; in particular, the cavity resonances are only slightly perturbed by 
the presence and position of the micromirror. The opposite is true of the |£| ^> 1 case, 
where there is coupling between pairs of cavity modes, typified by the avoided crossing 
in the spectra. The resonance frequencies can be obtained analytically, in the limit of a 
good bare cavity, as frequency shifts from the bare resonances: 
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with L c being the length of the cavity, x the position of the micromirror, and ko = 2n /\ 



the wavenumber of the light inside the cavity; Eq. (15) is identical to Eq. (4) in Ref. [25J. 



The two sets of solutions to Eq. (15) are, in the ( — > limit, separated by a free spectral 
range. These cavity resonances, plotted as detuned cavity lengths AL C = (L c /oj)Au, 
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Figure 3. Field intensity (left panels) at and friction coefficient (right panels) acting 
on the micromirror as the micromirror position (x) and cavity length (L c + AL C ) are 
scanned. The subfigures differ only in the polarisability of the mirror, as indicated. 
The cavity parameters are modelled from Ref. [12]. In the series of left panels, we note 
the progression from an almost bare cavity situation (a) to a very strong perturbation 
by the micromirror, leading to avoided crossings (f). The white dashed line traces a 
cavity node, whereas the black dashed lines [Eq. (JTs])] trace the cavity resonances. In 
the series of right panels, note that the friction coefficient is — as expected — a cooling 
force (blue) for red cavity detuning and a heating force (red) for blue detuning. The 
colourbars are on a logarithmic scale and are for 1 W of input power. 



are traced by means of the dashed black curves in the left panels of Fig. |3j We note 
that a unit on the vertical axis (AL C = A) is equal to twice the free-spectral range of 
the cavity. 

In the standard optomechanical coupling Hamiltonian, the mirror-field coupling is 
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Figure 4. Linear and quadratic optomechanical couplings as a function of mirror 
position for a very good cavity and for (a) ( = —1, and (b) ( = —10. In each 
figure we show the linear (solid curve) and quadratic (dashed curve) couplings, from 
Eqs. ( 18 ) and ( 19 1 . Note that the peak value of oj" is roughly proportional to £ whereas 
oj' is bounded. 



represented by a term of the form 

Hqm ~ f\hj'xa)a , (16) 

where x the position operator of the mirror, and oj' = d(Aoj)/dx. a is the annihilation 
operator of the field mode that has the dominant interaction with the micromirror; in 
the |C| — > limit, these field modes are the bare cavity modes of the whole cavity. 
However, as \(\ increases, the micromirror effectively splits the main cavity into two 
coupled cavities, giving rise to symmetric and antisymmetric modes, seen as the higher 
(bright) and lower (dark) branches in Fig. ^i) for < x < A/4; in such cases a is 
the annihilation operator belonging to one of these eigenmodes. We note that similar 
behaviour was observed in Ref. |25j . 

Certain effects, such as mechanical squeezing of the mirror position [32] and quantum 
non-demolition measurements on the mirror |33j, require not linear coupling to x but 
quadratic coupling to x 2 : 

i?Q 2 y ~ hJ'x 2 b)a , (17) 

with oj" = d 2 (Aoj)/dx 2 . In our notation, we have 

2k c ( sin(2k x) 
L c [1 + C 2 sin 2 (2%r)] ] 



and 



■ ^ r> u _\ 3/2 ^ i 



Ak^c ( cos(2/c x) 
[1 + C 2 sin 2 (2A; x)] ; 

One thing we note immediately is that there is no value for x such that oj' = oj" = 0; in 
other words, the optomechanical coupling is restricted to be linear or quadratic, to lowest 
order. Higher-order nonlinearities may be achieved by coupling different transverse 
modes of the cavity (see, e.g., the experimental results in Ref. [22]) but are overwhelmed 
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by the linear or quadratic couplings in a single-transverse-mode cavity. Moreover, the 
linear coupling u' is bounded in the ( — > oo limit: 

< ^ ~ 2tt x 8.42 MHz/nm , (20) 

with the numeric value corresponding to our parameters. In the same limit, u" exhibits 
resonant behaviour (see Fig. [4]), indicative of avoided crossings in the spectrum, peaking 
at a value of: 

->• ICI ~ 2tt x 0.10 |C| MHz/(nm) 2 . (21) 



L 



c 



We plot the lower (± — > — ) branches of Eqs. ( JX8| ) and (19) in Fig. [4] for two values for (: 
( = —1, representative of realistic micromirrors, and ( = —10, representative of a highly 
reflective micromirror. These correspond to cases (e) and (f) in Fig. [3j respectively. 
Coupling between the pairs of modes is not very strong for the ( = — 1 case; this is 
manifested by means of the smooth variation with x of co' and to" in Fig. |4]^a). The 
second case shows strong signs of the avoided crossing behaviour seen in Fig. [3]T), with 
oj' no longer behaving smoothly and oj" acquiring a resonance-like character. Note that, 
independently of the magnitude of (, the strongest quadratic coupling always occurs at 
the points where u' = 0. 

In parametrising our interaction in terms of a frequency shift Au we are effectively 
mapping the model originating from the TMM into a single-optical-mode model. It 
is important to note that this mode spans the entire cavity regardless of the nature 
of d what depends on ( is the spatial profile of the mode. In the limit ( — > 0, the 
field intensity is distributed uniformly throughout the cavity, whereas for large \(\, it 
is concentrated on one side of the membrane. These two situations are, as we have 
already discussed, handled differently in the CQED model, the former in terms of a 
single optical mode, and the latter in terms of two coupled optical modes. To highlight 
the failure of the coupled-optical-mode model as |£| decreases, we show in Fig. [5] the 
static force acting on the scatterer (i.e., the force when v = 0) as predicted by the two 
models. For the coupled-mode model, we use the predictions of Ref. (25], which hold 
for |r| — y 1, and deliberately misapply them to cases where |r| <C 1. From this model, 
given an input power P in , a tunnelling frequency g = c \t\ / L c , and a detuning A from 
resonance at x = 0, one obtains 

2lu'k c k 2 c + (A + Jxf - g 1 

he (2k c A) 2 + (k 2 c + cu' 2 x 2 + g 2 - A 2 ) 
with bJ = —2k c/L c . For large \(\, the two descriptions are essentially identical; indeed, 
it is easy to understand that the description of two coupled cavities is a good one when 
the reflectivity of the central mirror approaches or exceeds 90%. For reflectivities of the 
order of 50% (|£| ~ 1), however, the coupled-cavity description does not work well and 
one must switch to a scattering model to describe the situation accurately. For smaller 
Id still, as we have already mentioned, the predictions of the scattering model again 
agree with a CQED model of a scatterer (e.g., an atom) coupled to an unperturbed 
cavity field. 



p = ^T^-ru^ - y 
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Figure 5. Static force (i.e., the force acting on the mirror when v = 0) computed 
from the scattering model presented here (left panels) and a model based on a modal 
decomposition |25j (right panels), showing only one pair of modes. Red and green 
regions represent forces pointing in opposite directions, as indicated on the colourbar. 
We note qualitative agreement between the two models for x w 0.25A and for AL C close 
to the resonances, especially for large The discrepancies between the two sets of 
data, that are more pronounced for small polarisability, have significant consequences 
for any theory based on a coupled-cavity modal decomposition model. The black 



dashed lines [Eq. ( 15 1] trace the cavity resonances in the scattering model. The absolute 



values on the colourbar relate to the left panels. 
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4. Conclusion 

We have developed a generically-applicable theory to describe the motion of scatterers 
in electromagnetic fields. By applying this theory to the specific case of a scatterer in 
a cavity, we have shown how the scattering description can be used to bridge the gap 
between atom-CQED models, which rely on the atom interacting with one single mode 
that spans the entire cavity, and membrane-CQED models, where the membrane splits 
the cavity field into two coupled modes. It is in the region of current experimental 
interest, with membrane reflectivities of the order of 50%, that the discrepancy between 
the two descriptions starts emerging and where the usual "|r| — > 1" limit of membrane- 
CQED cannot be taken. 
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